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ABSTRACT 

We study the power spectrum of galaxies in redshift space, with third order perturba¬ 
tion theory to include corrections that are absent in linear theory. We assume a local 
bias for the galaxies: i.e. the galaxy density is sampled from some local function of the 
underlying mass distribution. We find that the effect of the nonlinear bias in real space 
is to introduce two new features: first, there is a contribution to the power which is 
constant with wavenumber, whose nature we reveal as essentially a shot-noise term. In 
principle this contribution can mask the primordial power spectrum, and could limit 
the accuracy with which the latter might be measured on very large scales. Secondly, 
the effect of second- and third-order bias is to modify the effective bias (defined as the 
square root of the ratio of galaxy power spectrum to matter power spectrum). The 
effective bias is almost scale-independent over a wide range of scales. These general 
conclusions also hold in redshift space. In addition, we have investigated the distor¬ 
tion of the power spectrum by peculiar velocities, which may be used to constrain the 
density of the Universe. We look at the quadrupole-to-monopole ratio, and find that 
higher-order terms can mimic linear theory bias, but the bias implied is neither the 
linear bias, nor the effective bias referred to above. We test the theory with biased 
N-body simulations, and find excellent agreement in both real and redshift space, pro¬ 
viding the local biasing is applied on a scale whose fractional r.m.s. density fluctuations 
are < 0.5. 

Key words: cosmology: theory - galaxies: clustering and redshift - galaxies: bias - 
large-scale structure of Universe 


1 INTRODUCTION 

The clustering properties of galaxies are dependent to a certain extent on the precise population analysed. This implies that 
not all populations can be unbiased tracers of the matter distribution - at least one must be biased. Simple bias models have 
been in existence since the early Cold Dark Matter (CDM) simulations suggested that pairwise velocities of galaxies would 
be too high unless the (optically-selected) galaxies were more strongly clustered than the mass (Davis et al. 1985). In the 
first simple models, it was assumed that the fractional overdensity in galaxies <5® = 5n/n was simply a multiple of the mass 
overdensity: 5® = b5. This linear bias model cannot be true in detail for all populations, since the shape of the power spectrum 
is unchanged in this case, and not all galaxy populations have the same shape of power spectrum, although the differences are 
not large (Peacock & Dodds 1994). At a more fundamental level, such a model could only survive if applied to a smoothed field, 
otherwise <5® < — 1 which corresponds to a negative galaxy density. The linear bias model can be viewed as an approximation 
to a more general Eulerian bias prescription, where the galaxy density is assumed to be some function of the present-day mass 
density (e.g. Coles 1993; Catelan et al. 1994; Weinberg 1995; Mann, Peacock & Heavens 1998). The linear bias term is then 
the first interesting term in a Taylor expansion of the function about 5 = 0. If the galaxy distribution is significantly biased 
with respect to the matter, there are major implications for cosmology, arising from the difficulty in detecting the matter 
distribution. On scales where masses can be reliably measured, such as in galaxy clusters, bias is required to reconcile the 
mass-to-light ratios with an Einstein-de Sitter Universe. Separately, linear studies of peculiar velocities, either directly (e.g. 
Sigad et al. 1998) or via redshift distortions (e.g. Hamilton 1992, Heavens & Taylor 1995; Cole, Fisher & Weinberg 1995; 
Hivon et al. 1995, Hamilton 1998, Hatton & Cole 1998) return the density parameter of the Universe only via /3 = Do ®/b, 
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so ignorance of b compromises determinations of ffo- In addition to the Eulerian bias discussed in this paper, there are other 
possibilities, such as Lagrangian bias (e.g. Davis et al. 1985, Bardeen et al. 1986) and stochastic bias (Pen 1997; Tegmark 
& Peebles 1998, Dekel & Lahav 1998), biasing determined by halo properties at some redshift (Mo & White 1996, Catelan 
et al. 1998, Catelan, Matarrese & Porciani 1998), or biasing determined by coherent processes over a large scale (Babul & 
White 1991; Bower et al. 1993). Despite a major industry in modelling (e.g. Kauffman et al. 1998), it is probably fair to say 
that we still have a great deal of uncertainty in precisely where galaxies should form, so it is not clear which, if any, of the 
above descriptions corresponds closely to reality. In this paper, we assume an Eulerian description for the bias, but we do not 
restrict it to linear form. 

The current generation of galaxy redshift surveys, such as the Anglo-Australian Telescope 2dF survey (Colless 1996), and 
the Sloan Digital Sky Survey (Gunn & Weinberg 1995) have sufficient space density and volume to look beyond linear theory, 
and this opens up possibilities of lifting the degeneracy between flo and 6, by measuring b empirically through higher-order 
statistics (Matarrese, Verde & Heavens 1997; Verde et al. 1998). Higher-order studies of the power spectrum have been made in 
real space by Makino, Sasaki & Suto (1992), Jain & Bertschinger (1994) and Baugh & Efstathiou (1994), and also in an elegant 
treatment in real and redshift space in the Zel’dovich approximation by Taylor & Hamilton (1996). Our treatment differs from 
previous perturbative calculations by its inclusion of bias and redshift distortions, and from the Zel’dovich approximation by 
inclusion of nonlinear, local bias and by the different treatment of nonlinear evolution. The layout of the paper is as follows: 
in Section 2 we set out the calculations formally, treating carefully the transition from real to redshift space; in Sections 3 and 
4 we separate the major effects which come into operation beyond linear theory, in real and redshift space, and in Section 5 
we present our conclusions. 


2 METHOD 


2.1 Real to redshift space 

Let r and s be the real and redshift-space coordinates, with the observer at the origin. The latter is defined to be the recession 
velocity (including the peculiar velocities v of galaxy and observer) from the observer, divided by the Hubble constant Hq. 
Ho = 1 is equivalent to using recession velocity as the distance coordinate, and we assume this from now on. Further we 
define pr(r) and ps(s) to be the density helds in real and redshift space. The mean density may be spatially-dependent, 
because of selection effects; we define the expected densities as <(>r(r) and <(>s(s). The overdensity in redshift space is dehned 
by 1 -I- 5s (s) = pa{s)/(j>s{s), and similarly for the real-space overdensity Sr- In all, we shall work with four random fields: 

• Sr'- mass overdensity in real space 

• Sa'- mass overdensity in redshift space 

• Sr', galaxy overdensity in real space 

• Sa : galaxy overdensity in real space 


and each will have its corresponding power spectrum Pf etc. For simplicity and consistency with previous notation, we drop 
the r subscript on the real space mass power spectrum. 

The coordinate transform from real to redshift space is (Kaiser 1987) 


s(r) = r 


^ ^ U(r) - 1/(0) 
r 


( 1 ) 


where U{r) = v • r/{Hor). Number conservation implies pr(i’)d^r = pa{s)(fis, which gives 

n , w M r ^ ^ [1 + ^’’(‘■)] 4>r{T) 

[1 + 54s)] (j)a{s) = i ^ 


where the Jacobian is 


J = 


1 + 


AU 


^ , dU{v)\ 


dr 


( 2 ) 

( 3 ) 


and AU{r) = U{r) — 1/(0). We make now the large-distant-volume approximation, where we assume that any modes we 
analyse have wavenumbers k which satisfy kr ^ 1 throughout. Terms AU/r, which are ~ S/{kr) if ~ 1, are ignored entirely 
in comparison with S. If we assume that <j) drops as some power of r, then a Taylor expansion of (f>{s) yields 4>{r) plus negligible 
correction terms. Beyond linear theory, the lowest-order contributions to the power spectrum arise from both second- and 
third-order terms in 5, so we expand to third order: 


Sa{s) = Sr{r) - U'{v) + U'\v) - Sr(v)U'{v) + Sr{r)U'\r) - U'Ar) = E(r), 


( 4 ) 


and ' = djdr. To linear order, differences in the argument do not matter, but to third order, they do. To the required order. 
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F(r) = F(s) - AUis)F'{s) [l - ;7'(s)] + iA;7'(s)F"(s), (5) 

which gives us the final expression for the redshift-space overdensity in terms of the real-space overdensity, all evaluated at s: 

Ss = Sr-U' + U'^ - SrU' + SrU'^ - U'^ - AUS', -f AUU" - SAUU'U" + 2AUU'5t + AUU''5r + ^AU^S” - ^AU^U”'. (6) 


2.2 Bias and evolution 


We assume that the galaxy overdensity field <5® = injn, where n is the number density of galaxies, is related to the mass 
overdensity field 5 via a local function. This plausible ‘local Eulerian’ bias model is chosen for simplicity and tractability, and 
has some support from simulation (e.g. Cen & Ostriker 1992, Kauffmann, Nusser & Steinmetz 1997). Of course other schemes 
are possible, such as Lagrangian bias (e.g. Catelan et al. 1998; see also Cole & Kaiser 1989, Mo & White 1996), or nonlocal 
bias (e.g. Bahcall & West 1992, Bower et al. 1993, Matsubara 1995). We expand the local function as a Taylor series around 
5r = 0 (Fry & Gaztanaga 1993): 




1=0 


r- 


(7) 


An unbiased galaxy field would have 6i = 1 and all other bias parameters zero. Already there is a subtlety. We will truncate 
the expansion, which will only be a good approximation if the value of 5 is typically much less than unity. Therefore our 
biasing assumption is that there is some smoothing scale for which (j^ is a good approximation, and the fields above should be 
interpreted as the smoothed galaxy and mass density fields. For the perturbative expansion to be valid, we choose a smoothing 
scale large enough that the smoothed density field has small fluctuations (this will be quantified later). If we wish to smooth 
on a smaller scale, then a numerical approach is probably necessary (Mann, Peacock & Heavens 1998). Note also that the 
galaxy distribution is of course a point process. We assume that the positions are determined by a Poisson sampling of a 
density field whose overdensity is &g. We do not consider stochastic biasing (Pen 1997), and ignore shot noise in the power 
spectrum. 

The next non-zero terms beyond linear theory will be of order 4 in 5r, so we need to keep terms up to j = 3. Throughout 
we ignore the bo term, required to ensure that (5®) = 0, since we will be interested in the spectral properties of the galaxy 
field, and drop all terms which contribute only to k = 0. 

Inserting (|^ into (^), we get, to third order in Sr’. 


<5? = [biSr - U'] + 


y (5^ -b t/'^ - biSr U' - biS'^AU -b AU U'‘ 


+ 


- -j-slu' + blSr U'^ - U'^ - b2AUSrS'r 


-3AU U'U" + 2biStAU U' + b^SrAU U" + ^biS'^{AUf - ^{AUfU" 


( 8 ) 


The Fourier transform of this may be taken in redshift space, making use of the transforms of simple products: 

(A:y)k = J TiiiTk2S°{k-ki--k2)XiY2 

(Xy^)k = / d®kld®k2d®k35^(k-kl-k2-k3)Xiy2^3 

( 27 r)'> J 

(9) 

where Xi is the ki component of Ak etc, and 5^ is the Dirac delta function. In the distant-observer approximation, 8/dr —> 
ikfi, where /r = k • r/k and r is a unit vector from the observer to the galaxy (assumed constant across the sample) and the 

transform of U is Uu = where —r^k is the transform of the velocity divergence. / = d\nD/d\na ~ Do'®i where D[a) 

is the growing-mode amplitude and a is the scale factor of the Universe. Second-order terms were computed by Peebles (1980); 
the third-order expansion of the fluid equations (see Fry 1984) is detailed in Catelan & Moscardini (1994a) and Catelan & 
Moscardini (1994b), giving 5k: 

5k = ek-b^^y'd"kid®k25^(k-ki-k2)4"^ki,k2)eie2 

+ /d®kid3k2d®k35°(k-ki-k2-k3)J®(ki,k2,k3)£ie2e3. 

(27r)*> J 

( 10 ) 

ck is the linear, real-space Fourier coefficient of the density field 5(r). The corresponding expressions for r/u, involve replacing 
J by K. The functions J and K are quoted for Dq = 1, as computed by Goroff et al. (1986), Gatelan & Moscardini (1994a) 
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and Catelan & Moscardini (1994b): 

d 2 )n 1 ^ 5 , ki • k 2 /fci fc 2 ^ , 2 ^ki • 

/fci k2\ 4 /ki • k 2 \^ 

U 2 fcJ 7 V fcife ) ’ 


(ki,k 2 ) 


5 ki ■ k 2 

7 2fcifc2 
3 ki ■ k 2 

7 2kik2 


J®(kl,k2,k3) 


_ 7 ( 2 ) 


J^"'(k2,k3) 


1 1 ki • (k 2 + k3) 4 k • ki k ■ (k 2 + k 3 ) 
3 3 (k 2 + k3)2 9“fcf (k 2 + k3)2 


2 k • ki k • (k 2 + k 3 ) ks • (k 2 + k 3 ) 1 k • k 2 k • ks 


9 kl (k2+k3)2 


fc| 


9 fci 




7^®(kl,k2,k3) = 3/"^(kl,k2,k3) - j^J<"^(k2,k3) - j^<")(kl,k2) 


kf 


(kl + k 2)2 


( 11 ) 


with k = kl + k 2 + k 3 in the last two expressions. The subscript S indicates that the expression has been made symmetric 
w.r.t. any permutation of the arguments. If not, then the symmetrized kernel must be obtained by averaging the quoted 
expression over all permutations. These kernels are correct for flo = 1 but are only weakly dependent on Qo (e.g. Bouchet 
et al. 1992, Catelan et al. 1995; see also Bernardeau 1994b, Eisenstein 1997, Kamionkowski & Buchalter 1998). After some 
manipulations, the transform of the third-order, biased, redshift-space density field is 




— pC) 


FW(k)ek + ^^ J d^kid®k25^(k-ki-k2)E|">(ki,k2)eie2 

f d®kid®k2d®k35^(k - kl - k2 - k3)E®(ki,k2,k3)eie2e3 


( 2 . 


with the following kernels: 


( 12 ) 


= 


Ef(kl,k 2 ) = 


bi + fv 

61 jf (kl, k 2 ) + (kl, k 2 ) + ife + ^ 


2 , 2 , 1 

fki k2V 

Ml + P 2 + M1M2 1 

U 2 kl). 


+f 


2 2 , 

niV2 + -TT- 


^ 2 kl 2 ^2 ' 


E®(kl,k2,k3) = 


&iJ^®'(ki,k2,k3) + /p^A'^®^(ki,k2,k3) -I- + ^ + 


r2 .2^1 , L J.2,. .3^3 , 61j.2..2 V 


+ 2bif ^iip2/i37r + bif P2M37r + T-C , , 

K 2 K 2 2 K 2 K 3 

+ /"^(k2,k; 


1 I /•3 2 2 2 I n /•3 2 3 ^3 . 1 _f3 d V 


, ro ^ ^ ^ , n SO ^ o •'"O , ro 4 

+ / MiA^2M3+3/ + 2 / 


) (b2 + biffli + fei/^ii/i2+3^|^) 

A'*'^^(k2, k3) ^5l//i2+3 + blffll)J,2+3 + 2/^/if/ii+3 -I- 12 . 1 ^ 2 + 3 —+ f 


2 3 kl 

^1^2+3- 


k2+3 


(13) 


A subscript S on the F terms also indicates that the term has been symmetrized w.r.t. its arguments. Note that the last term 
has not been symmetrized and jj,2+3 = (k .2 -b k 3 ).f/|k 2 + ks]. 

From these equations we can obtain the power spectrum to third-order , which includes non-zero correction terms of two 
types, as in the real-space unbiased case (Makino, Sasaki & Suto 1992; Jain & Bertschinger 1994). The redshift-space power 
spectrum (k) is defined by 

= (27r)^Pi'(ki)J°(ki + k2) (14) 

and the real-space mass linear power spectrum Pii{k) is dehned similarly by (ekiek 2 ) = (27r)®Pii(fci)J^(ki + k 2 ). The 
Gaussian nature of the initial fluctuations implies, by Wick’s theorem, that (£ 16263 ) = 0 and (£i£ 2 £ 3 £ 4 ) = (£i£ 2 )( 6364 ) plus 
cyclic permutations. Hence we hnd two (one-loop) terms which arise at the next level beyond linear theory for the spectrum 
(tree level): 

Pi'(k) = Pi'll+Pf22+n®13 


= {l+P^^^fblPll{k)+2 


d®q 

(27r)3 


Pii(?)Pii(|k-q|) 


Pf (q,k-q) 


+ 6(1-b)6iPii(fc) 


(Vq 


(271 


Pii(<7)P® (q,-q, k). 
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(15) 

where f3 = // 6 i. In the unbiased case in real space {f3 = 0), these formulae reduce to those in Jain & Bertschinger (1994) 
(note there is a factor of missing in their equation 19). For this case, the limits of integration for the P 22 ^13 terms are 

not too important, provided the power spectrum decays sufficiently rapidly at high k (faster than k~^), and this is fine for 
spectra of practical interest. For the bias terms, however, one has to be careful, because if one places a cutoff in the integrals, 
one finds that the resulting power spectrum is sensitive to this cutoff. This should not be surprising: the cutoff effectively 
smoothes the linear power spectrum (with a top-hat filter in fc—space), and the power spectrum depends on what scale the 
filtering takes place. Furthermore, it makes physical sense: if a local bias prescription applies at all, it will be necessary to 
specify on what scale it acts. These considerations lead us to the following Eulerian bias proposal: 

• Evolve 5 to 3rd-order 

• Smooth with some filter function corresponding to a smoothing scale Rf 

• Apply the Taylor expansion to third-order to get Sg. 

The resulting field is evidently itself smoothed. This procedure is effected by replacing all occurrences of the linear by e^Wk, 
where Wk is the transform of the filter function [e.g. Wk = exp(—A:^i?^/2) for a Gaussian filter]. We shall not explicitly retain 
these factors of Wk', the smoothing is implicit in tk- Note that the operations of smoothing and biasing do not commute; 
unless the matter power spectrum has a natural small-scale cutoff, the order assumed here is the only sensible choice for a 
perturbative calculation. This is a critical assumption, of course, but we reiterate that it has some support from simulation 
(e.g. Cen & ©striker 1992). Eor further discussions of the effects of smoothing, see Coles, Melott & Shandarin (1993), Melott 
& Shandarin (1993), Bernardeau (1994a). 

2.3 Bispectrum 

The galaxy redshift-space bispectrum Bf(ki,k 2 ,k 3 ) is the 3-point function in Fourier space, defined by 

(<Iski<Isk 2 'Isk 3 ) = (27r)®Bf (ki, k 2 , k3)(5^(ki + k 2 -b k 3 ). (16) 

In real-space, this has been studied by Matarrese, Verde & Heavens (1997), Scoccimarro (1997) and Scoccimarro et al. (1998). 
We can use the preceding formalism to obtain the biased, redshift-space bispectrum, in the distant observer approximation: 

Bf(ki,k 2 ,k 3 ) = 2 Pii(fci)Pii(fc 2 )EW(ki)E«(k 2 )E^"^(ki,k 2 ) -b cyclic terms (1,2) ^ (2,3) and (3,1) 

= 2Pn(fci)Pn(fc2)&?(l + /I/r?)(l -b/lM2) (ki, k2) -b/ 3 /i 2 A( 2 '(ki, k2) + ^ + bil3^+ 

L Zui 2 

biP / ki k2\ bl 2 f 2ki 2 fc 2 \l 
-b —biiM 2 -b — J -b —f3 fiiU 2 + U 2 -j^ j + cyclic terms. 

(17) 

In a separate paper (Verde et al. 1998) we use this to investigate how redshift surveys can be used practically to estimate 61 
and hence determine the density parameter (via /3). 

3 NEW FEATURES: REAL SPACE 

Interestingly, there are some new features which arise in real-space from biasing alone. Consider the real-space corrections to 
the power spectrum, Pfj = Pa 22 {P ~ 0) 8 - 11 ^ P 13 = Pfsih = 0). At small fc, the former is more-or-less constant with k, while 
the latter is proportional to Pii(fc). We look at each of these in turn. 

3.1 Constant power on very large scales 

On very large scales, where A^(fc) <C 1 (where is the contribution to the variance in mass overdensity per unit In(fc)), the 
terms arising from nonlinear evolution (and a linear bias 61 ) are generally small (Jain & Bertschinger 1994), so we concentrate 
on the contributions from nonlinear bias. 

The P 22 term is 

^ 2 ® 2 (fc) = ^ /d"klP(kl)P(|k-kl|) [ 6 ljf)(kl,k-kl) + y]" 

(18) 

which provides a constant contribution to the power when |k| is small. The presence of such a term was suggested by Coles 
(1993) and Strauss & Willick (1995); here we provide a mechanism for calculating it perturbatively for arbitrary local bias. 
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k 

Figure 1. Real-space power spectrum. The (upper) solid line is the final power spectrum; the dashed lines are the linear power spectrum, 
unsmoothed and smoothed, and the dotted lines are (from the bottom at fe = 10“®) the contributions to the power from P 22 , P 13 and 
their sum. The dot-dashed line is the second-order mass power spectrum, which merges with the dashed linear power spectrum to form 
the lower solid line. For model details, see text. 


For an underlying matter power spectrum which falls to zero as fe ^ 0, P 22 will eventually dominate over both the 
linear term and P^g, leading to a divergent effective bias = P®(fc)/P(fc) in this limit (P is the matter power spectrum 

evolved to third order; on these large scales, linear theory is accurate, so P ~ Pii). In turns out, for an underlying Zel’dovich 
spectrum, this term is unlikely to be of importance unless the other correction term below is small. The second-order 
power spectrum is shown in Fig. 1 for a CDM-like spectrum with shape parameter F = 0.25, with and the amplitude is chosen 
such that A^(fc = 0.1) = 0.03. The mass field is smoothed on a scale of 2 h~^ Mpc, so that the linear smoothed variance is 
o-Q = 0.20, and the bias parameters are 61 = 62 = 63 = 1- Note that the nonlinear bias masks the effect seen in the unbiased 
case by Jain & Bertschinger (1994) of a transfer of power from large to small scales. 

The constant power at large scales should not be confused with the requirement that P(fc) ^ 0 as fc ^ 0, which applies if 
the mean density is estimated from the survey. This will suppress the power for k <1IL, where L is the characteristic depth 
of the survey, and is quite independent of the considerations here. Similarly, the neglect of constant terms in <5 affects only 
k = 0. 

The term here has been conjectured before, and it is appropriate to make a few remarks about its nature. It is essentially 
a shot-noise term arising from the peaks and troughs in the underlying density field being nonlinearly-biased by the quadratic 
term 62 • The biasing gives a contribution to the density field consisting of spikes at the peaks and troughs. We consider here 
the case when the underlying matter power spectrum tends to zero on large scales, so we look only at this remaining ‘spike’ 
contribution (see Appendix A). It may seem odd that the troughs make a contribution; is should be remembered that they 
do not necessarily correspond to actual overdensities, as there is a linear contribution which generally lowers the density - 
a positive quadratic term merely means the underdense regions are not quite as underdense as a linear bias would imply. 
However, by assumption, this linear contribution has zero power on very large scales, so the quadratic term dominates. 

To illustrate the shot-noise nature of the spikes, consider a non-pathological power spectrum, filtered to remove power 
with k > fcmax- It is straightforward to show that the constant-power term is roughly 

P^2-blA\k,„,„)- (19) 

n 

where A(fcmax) measures the nonlinearity on the smoothing scale, and n ~ fc,0ax is the number density of peaks and troughs on 
this scale (Peacock & Heavens 1985, Bardeen et al. 1986). Shot noise of a point process gives Pshot = 1/h (e.g. Peebles 1980); 
the prefactor arises here because each spike has a height of roughly & 20 'o so the integral under each spike is ~ 62 A^(feniax)- 

This shot noise argument would be dangerous if it was applied to peaks alone, as we know that these are strongly 
clustered (Kaiser 1984, Peacock & Heavens 1985, Bardeen et al. 1986), and their clustering might well exceed the shot noise 
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Figure 2. Effective bias in real space, from perturbation theory (solid). This is approximately constant over a fairly wide range of k, 
yet the bias is far from linear. The dotted line shows the approximate analytic formula (^). The filter erases power at k >0.2 and the 
spike at fc = 1 is due to the matter power spectrum crossing zero, by which time third-order perturbation theory has broken down. 


contribution. In the appendix we prove that high peaks and low troughs, taken together, are very weakly clustered, so the 
power contribution is indeed essentially a shot noise term from unclustered small-scale spikes in the biased field. 


3.2 Effective bias on large scales is not bi 

If the P 22 term can be neglected in comparison with Pfg, then the effect of biasing and second-order evolution is to add a 
term which is proportional to the linear P. This holds over a wide range of scales, but not at the very largest scales, where 
the power spectrum is dominated by the constant contribution from Pj'j, leading to a divergent bias. 

The fact that Pfg is proportional to Pn gives a constant effective bias = P®(fc)/P(fc) ~ P®(fc)/Pii(fc)]. In linear 
theory, the effective bias is simply 61 ; in this section we compute the corrections to this value. In real space, only 3 terms 
of P^®^ survive. The first, proportional to comes from nonlinear evolution of the mass field, and is very small on large 
scales, and we neglect it here. There is a straightforward contribution from fea/fi, which gives 

bib3Pii{k)aQ ( 20 ) 


where ctq = (27r^) ^ J dkk^ Pii(fe) is the linear variance in the smoothed field. The other nonzero term is 


2b2 

(27r)= 


d®kPii(q) 


Js^ (-q, k) -b (q, k) + Tf’ (q, -q) 


7 ( 2 ) 


( 2 )/ 


( 21 ) 


The last term in brackets is zero because J is zero, and the first terms are equal. Integration over angles gives a contribution 
again proportional to ctq, and the final result is 

Pf 3 (fc) = felPll(fc) + 63) (22) 

The Pfg term is of more practical importance than Pjj. We see that it is proportional to the linear power spectrum, so it 
corresponds to a scale-independent bias, which, when including the linear terms, gives an effective bias of 

b ! ff ( k ) = 6? + 61 + 63) <^1 (23) 

The level of bias depends on the variance of the smoothed field, for which the truncated Taylor expansion is assumed to be 
a good description of the galaxy density field. It also depends on both the second- and the third-order bias coefficients. The 
suggestion that the bias should be constant on large scales has been made by a number of authors (Coles 1993; Scherrer 
& Weinberg 1997; Mann, Peacock & Heavens 1998), and we confirm that this is true to third order in perturbation theory, 
provided one is not looking at the very largest scales where P^j dominates. 

The effective bias is plotted in Fig. 2, along with the prediction of (^^, which is seen to be a good approximation over a 
wide range of fc, for these parameters. In Fig. 3 we demonstrate that the approximate formula works well for a biased N-body 
simulation, with the same linear power spectrum, normalised to ag = 0.26, from the Hydra consortium (Couchman, Thomas 
& Pearce 1995). The rather small amplitude of the power spectrum gives us a large dynamic range between the smoothing 
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Figure 3. The effective bias factor Pbiased!Pmass for a Hydra N-body simulation (solid), along with the approximate formula 

for bi = b 2 = bs = 1 and Rf = 2h~^ Mpc, which gives an r.m.s. fractional overdensity of 0.41. For other details, see text. This figure 
also differs from Fig. 2 in that the biased field has been deconvolved with a Gaussian, which affects the results at high k. 

scale and the size of the simulation. The smoothing is done with a Gaussian filter of radius 1 h~^ Mpc, giving a variance of 
(Jq = 0.17, and we take &i = &2 = 63 = 1. Note that a bias function which includes non-zero bi terms beyond i = 3, such as the 
exponential high-peak biasing exp(a5) (Politzer & Wise 1984, Jensen & Szalay 1986, Bardeen et al. 1986) gives essentially 
the same result provided that the variance in the field is not too large. For the power spectra here, this requires aao < 0.3, 
but the exact figure depends on how skewed the field is. As the field becomes progressively more nonlinear, the bias increases 
slightly towards high k, and the approximate formula begins to underestimate the bias. Note that the simulations do not 
probe a large enough scale to detect the constant contribution to the power. 

Since Pfg is proportional to the linear power, it is formally possible at this level of perturbation theory to cancel the 
linear power altogether by appropriate choice of 62 and 63, leaving the constant P 22 term. Perhaps not surprisingly, it appears 
impossible to do this on N-body simulations without perturbation theory breaking down. 


4 REDSHIFT-SPACE 

In Fig. 4 we show the effects in redshift space of biasing and evolution, for different angles of wavevectors with the line-of- 
sight, showing similar qualitative behaviour to Kaiser’s (1984) analysis. In Fig. 5 we show the redshift-space power spectrum 
(averaged over angle) for the Hydra simulation, for the same smoothings and parameters as in Fig. 3. We see here remarkably 
good agreement with perturbation theory of biased fields in redshift space. Note that linear theory, with linear bias, predicts 
that the redshift-space and real-space power spectra should be proportional to each other. Here we see that the higher-order 
terms give very different behaviour. Note also that perturbation theory gives excellent results all the way down to Mpc scales, 
where the smoothing is applied. 


4.1 Quadrupole-to-monopole ratio 


Hamilton (1992) proposed the quadrupole-to-monopole ratio of the power spectrum as a diagnostic for /3. These are obtained 
by expansion of the power spectrum in a shell of fixed |k|, in terms of Legendre polynomials (in p) of order 0 and 2. In linear 
theory, the ratio of the coefficients "P? has the value 


n 


jP + 

1 + 1 / 3 +i/32 


(24) 


Taylor & Hamilton (1996) extended this idea by using the Zel’dovich approximation to investigate the behaviour of the 
quadrupole-to-monopole ratio into the mildly nonlinear regime. In this section we use perturbation theory results to investigate 
the influence that bias might have on conclusions drawn from the ratio. In practice, it is best to use the full 3D power spectrum 
to estimate parameters, but analysis of the quadrupole-to-monopole ratio can be a useful aid to see why we might be able to 
extract information, and why there might be problems. In Fig. ^ we see that, on large scales, linear theory accounts well for 
the ratio, with significant deviations appearing on small scales in qualitative agreement with Taylor & Hamilton (1996). The 
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CO 



Figure 4. Redshift-space power spectrum for the model of Fig. 1. The curves are for values of p from 0 (bottom) to 1 in steps of 0.1. 
The dashed and dotted lines are the unsmoothed linear and nonlinear power spectrum in real space. 





Figure 5. Hydra simulations in redshift space (points), along with the theoretical curve from perturbation theory (fei = 62 = 63 = 1.0, 
and (70 = 0.41). The dotted line is the smoothed linear theory matter power spectrum in real space. 


slight discrepancy with linear theory may arise from only computing the power spectrum with a rather crude separation of 
A/r = 0.1. On small scales, the qnadrnpole-to-monopole ratio was studied by Taylor & Hamilton (1996) and Hatton & Cole 
(1998), where it is affected serionsly by caustic formation and small-scale virialised structures. We confine our remarks to 
larger scales (k <0.1), where these effects are not important, but there are still some surprises. Fig. shows the quadrupole- 
to-monopole ratio for a biased CDM-like held, with F = 0.25, with 61 = 62 = 63 = 1-0 and ctq = 0.196 and Rf = 2h~^ Mpc. 
There is a large range of wavenumbers for which the ratio is very nearly constant, which might give one confidence that a 
linear bias is an adequate description. However, the retrieved value of /3 — 0.81 agrees neither with linear theory (/3 = 1), 
nor with the appropriate from the effective bias on large scales (/3 = 0.74). It is the latter which we would like to estimate 
accurately, since it is the effective bias which tells us the amplitude of the underlying matter power spectrum. In this case, 
the recovered effective bias (1.24) is not close to the true effective bias (1.35). The discrepancy arises because the one-loop 
corrections do not have the same angular dependence as the Kaiser factor multiplying the linear power. This is not really 
a breakdown of linear theory as such; the underlying mass field is still well described by linear perturbation theory. Here it 
is the nonlinear bias which causes the problem, and the simple analysis of Kaiser (1987) cannot be applied. These findings 
mean that attempts to measure /3 from linear redshift distortions (e.g. Hamilton 1992; Fisher, Scharf & Lahav 1994; Heavens 
& Taylor 1995; Ballinger, Heavens & Taylor 1995; Cole, Fisher & Weinberg 1995; Tadros et al. 1998) must make the further 
assumption that the one-loop corrections are small. 
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in 



k 


Figure 6. Quadrupole-to-monopole ratio for the mass, as a function of wavenumber (solid). Beyond k = 4.5, the filtered power is zero 
to machine accuracy. The dashed line is the linear theory ratio (M). 


in 



k 

Figure 7. Quadrupole-to-monopole ratio of the biased field (solid), as a function of wavenumber, along with unbiased linear theory 
(dotted) and linear theory, but with b replaced by fegfj from perturbation theory (dashed). The behaviour at fc < 10“® arises from the 
constant-power term in the biased field. 


5 CONCLUSIONS 

In this paper, we have shown how the power spectrum is altered in real and redshift space by the effects of nonlinear local 
bias and evolution. The assumptions in the analysis are that the smoothed galaxy density field is some local function of the 
underlying smoothed present-day mass field (i.e. Eulerian local bias). The smoothing scale is assumed to be large enough 
that a Taylor expansion to third order: 5g ~ X/j-o suffices. The approach we have used is perturbative; the one-loop 

corrections give contributions from third-order (in ^) which are as important as second order, so we have kept evolution 
and bias terms to this order {j = 3). We confirm the predictions of Scherrer & Weinberg (1997) that the power spectrum 
has a constant component, plus a roughly scale-independent bias on large scales. We show that the constant power term is 
essentially a shot noise term arising from quadratic biasing of peaks and troughs. The most important results in real space 
are the expressions for the magnitudes of the two effects: equation (|^) for the effective bias on large scales, and (^|) for 
the constant contribution, which could in principle make difficult the task of unveiling the underlying power spectrum on 
extremely large scales. However, our estimates of the size of this effect suggest that it is not going to be a practical problem 
in the foreseeable future, since it is likely to dominate only on scales of the order of the horizon size. Should the galaxy power 
spectrum be measured on extremely large scales at some time in the future, there is a possibility that a deviation from the 
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scale-invariant spectrum (in |fc|) might be interpreted as evidence for an open Universe (Kamionkowski & Spergel 1994). This 
paper provides an alternative interpretation in terms of nonlinear bias. 

Potentially more serious is the effect of nonlinear bias on the amplitude of the power spectrum, in both real and redshift 
space. Here there is the possibility of confusion, in that the shape of the linear power spectrum is preserved over a wide range 
of k, so one can define a scale-independent effective bias from the ratio of the galaxy to matter power spectra. The difficulty 
is that this bias is not the linear bias (6i in the Taylor expansion of the galaxy density held), nor is it the bias which would be 
recovered from studying anisotropy in the redshift power spectrum. These complications may compromise efforts to deduce 
cosmological parameters from redshift distortions or peculiar velocities, unless a further assumption is made that the one-loop 
corrections are small. Without a better understanding of galaxy formation, it is impossible to make this assumption with a 
great deal of conhdence. What would be interesting to see is whether current models can identify a scale on which an Eulerian 
bias is a good description of the galaxy held; if such a scale can be identihed, then the computations in this paper may be 
useful in providing analytic approximations for the effects on the power spectrum and redshift distortion. 
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APPENDIX A: CORRELATIONS OF HIGH AND LOW REGIONS 

A series expansion for the correlation function of regions above a threshold was obtained by Jensen & Szalay (1986); see also 
Poiitzer & Wise (1984). The first term in this series is the one obtained by Kaiser (1984), which shows that high regions 
can be very strongly clustered. In this appendix, we calculate the correlation function of high regions (above a threshold) 
plus troughs (below a threshold), and show that high regions plus troughs are very weakly correlated on large scales, with a 
correlation function which is proportional to the square of the correlation function of the field. 

We follow the formalism and notation of Jensen & Szalay (1986). Let the variance in the Gaussian field be cr^, and 
consider high regions above S = va and troughs below 5 = —va. 

The two-point correlation function fx of regions separated by r, with v in some region X is given by an integral over the 
normal bivariate Gaussian: 

^2;^ 

where 

1= f -^exp(-:/V2) (A2) 

JX V JtT 

and the covariance matrix Cij = {viVj) is unity on the diagonal, and ip = the diagonal. 

To evaluate the double integral, we use the same trick as Jensen & Szalay (1986). We take a Fourier transform of the 
exponential: 


1 ^-1 


(Al) 
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r 1 / 1 \i 

f d^k , 

■ 1 

exp 

-2 

= / --— exp(— ik • v) exp 

/ ZTT 

— 2 (kiCijkj) 


Then the double integral may be written 

f f f (fk 

J J dvidv 2 J —exp(—ik • exp(—fcifci/2) exp(Q), 

where Q = —kik 2 ip. We expand exp(Q) = lm\, and then regard Q as an operator 

d d 

Qexpf— ik • v) = lb --——— exp(— ik • v). 

OUl dV2 

The double integral can then be written 
j j du,du 2 j^^(£-J-^ exp [-(.? + .|)/2]. 

X J X n 


(A3) 


(A4) 


(A5) 


(A6) 


Jensen & Szalay (1986) took the region to be Vi > v. We take the region to be |t'i| > v, which gives four contributions to the 
double integral, and 7 is also modified to twice the high region result, i.e. erfc(i^/\/2)- The terms in the double integral from 
peak-peak and trough-trough give the same as Jensen & Szalay (1986), i.e. 


2 ^ ^ [77,„-i(i^)exp(-iz^/2)]^ , (A7) 

m=0 

where the Hermite polynomial is H^ix) = exp(a;^/2)(—d/da;)'” exp(—a;^/2). The cross-terms in the double integral give an 
identical contribution, but multiplied by (—I)*”. The effect of this is to eliminate the odd terms in the series, whilst leaving 
the even terms intact. The strong clustering term (for high regions alone) which Kaiser computed disappears, leaving the first 
non-zero contribution as 

Cx = y + 0(^") (A8) 

so if the underlying density field has small correlations on large scales (-!/> ^ 0), the field consisting of high regions and low 
troughs is very weakly clustered. 
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